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ABSTRACT  OF  RESEARCH  ACCOMPLISHED 
ON  CONTRACT  N00014-79-C-0522 


The  research  performed  on  this  contract  addressed  statistical 
methodologies  that  would  be  appropriate  in  censored  survival  analyses 
and  competing  risks  models.  Both  general  parametric  and  nonparametric 
models  were  considered.  The  problems  of  efficient  estimation  of 
parameters  of  the  underlying  failure  time  distribution  were  addressed 
when  the  data  may  only  be  partly  observable  due  to  the  presence  of 
censorship.  The  methodogy  developed  permits  costs  associated  with 
monitoring  of  experimentation  and  of  recruitment  of  units  to  be 
incorporated  into  the  basic  infrastructure  leading  to  the  description 
and  properties  of  optimal  procedures  that  minimize  overall  expected 
cost.  These  optimal  designs  are  inherently  sequential  in  nature,  in 
that  the  optimal  stopping  time  or  sample  size  is  solely  determined 
by  the  observed  information  gathered  by  the  experimeter. 
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FINAL  REPORT 
on 

ONR-Contract  N00014-79-C-0522 
by 

Joseph  C.  Gardiner 
Principal  Investigator 

This  report  describes  the  main  achievements  of  our  research  efforts 
on  the  above  contract.  The  central  thesis  of  our  research  related  to 
statistical  inference  problems  arising  in  the  analysis  of  censored  data. 
Such  data  or  incomplete  observations  arise  naturally  or  by  design  in 
several  longitudinal  studies  which  encompass  clinical  trials  and  life¬ 
testing,  reliability,  demography  and  epidemiological  investigations.  Our 
research  publications  cover  three  broad  areas:  i)  Sequential  estimation 
procedures  from  randomly  censored  data  ii)  Properties  of  likelihood 
ratio  processes  arising  in  progressively  censored  data,  and  iii) 
Estimation  of  biometric  functionals  in  multiple  decrement  models. 

Listed  below  are  the  publications  prepared  under  this  contract. 

It  is  followed  by  a  brief  description  of  their  contents  and  significance 
of  the  results  achieved. 

a)  List  of  publications 

1.  Risk-efficient  estimation  of  the  mean  exponential  survival  time 
under  random  censoring.  Proc.  Natl.  Acad.  Sci.  USA,  Vol.  81,  pp.  5906- 
5909  (1984).  J.  Gardiner  and  V.  Susarla. 

2.  Sequential  estimation  of  the  mean  survival  time  under  random 
censorship.  Common,  Stat.  Sequential  Analysis  Vol.  2,  201-233  (1983). 

J.  C.  Gardiner  and  V.  Susarla. 

3.  Fixed-width  confidence  intervals  of  the  exponential  mean  survival 
time  under  random  withdrawals.  J.  Stat.  Plann.  &  Inference.  12.  153-159 
(1985).  J.  C.  Gardiner,  V.  Susarla  and  J.  van  Ryzin. 
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4.  Time-sequential  estimation  of  the  exponential  mean  under  random 
withdrawals.  To  appear  in  Ann.  Statist.  (1986).  J.  C.  Gardiner,  V.  Susarla 
and  J.  van  Ryzin. 

5.  On  the  estimation  of  the  median  survival  time.  To  appear  in 

Inst.  Math.  Stat.  Mono.  &  Lecture  Notes  (1986).  J.  C.  Gardiner  &  V.  Susarla. 

6.  On  the  asymptotic  theory  of  estimation  of  the  survival  function 
in  competing  risks  analyses.  Michigan  State  University  Research  Mono. 

(1984).  J.  C.  Gardiner  and  V.  Mandrekar. 

7.  Weak  convergence  of  a  Bayesian  nonparametric  estimator  of  the 
survival  function  under  progressive  censoring.  Statistics  &  Decisions 
1,  257-268.  J.  C.  Gardiner  and  V.  Susarla. 

8.  A  nonparametric  estimator  of  the  survival  function  under  progressive 
censoring.  Inst.  Math.  Stat.  Mono.  &  Lecture  Notes.  Vol .  2,  26-40.  (In 
Survival  Analysis.  J.  Crowley,  R.  Johnson,  editors).  J.  C.  Gardiner 

and  V.  Susarla. 

9.  Local  asymptotic  normality  of  progressively  censored  likelihood 
ratio  statistics  and  applications.  J.  Multivariate  Analy.  12,  (2)  230-247 
(1982).  J.  C.  Gardiner. 

10.  Convergence  of  progressively  censored  likelihood  ratio  processes 
in  life-testing.  Sankhya  V43  I,  Part  A,  37-51  (1981).  J.  C.  Gardiner. 

11.  Asymptotic  normality  of  a  variance  estimator  of  a  linear  combina¬ 
tion  of  a  function  of  order  statistics.  Z.  Wahrsheinl ichkeitstheorie 

und  verw.  Gebiete  50,  205-221  (1979).  J.  C.  Gardiner  and  P.  K.  Sen. 

12.  Properties  of  time-sequential  statistics  in  life-testing. 

Reliability  &  the  Acquistions  Process,  1983-193.  D.  DePriest,  R.  Launer, 
editors.  J.  C.  Gardiner. 


13.  Asymptotic  normality  of  class  of  time-sequential  statistics 
and  applications.  Common.  Stat.  A7,  4  373-388  (1978).  J.  C.  Gardiner 
and  P.  K.  Sen. 

14.  The  asymptotic  distribution  of  mortality  rates  in  competing 
risks  analyses.  Scandinavian  J.  Stat.  Vol.  9,  Pt.  1,  31-36  (1982).  J. 

C.  Gardiner. 

15.  Properties  of  occurence-exposure  rates  in  multiple  decrement 
theory.  J.  Stat.  Plann.  &  Inference.  Vol.  8.  301-314.  J.  C.  Gardiner. 

16.  Bounded  length  sequential  confidence  intervals  for  the  median 
survival  time.  Submitted  for  publication. 

The  following  titles  are  in  preparation. 

17.  A  comparison  of  sequential  confidence  interval  estimation  procedures 
for  the  median  survival  time.  A  Simulation  Study. 

18.  Risk-Efficient  estimates  for  quantiles  of  the  survival  distri¬ 
bution  with  sampling  costs. 

19.  The  boundedness  of  the  regret  of  a  sequential  procedure  for 
the  estimation  of  the  exponential  mean  failure  time. 

20.  Sequential  point  estimation  of  the  mean  residual  life  under 
random  censorship. 

21.  Fixed-width  confidence  intervals  for  the  regression  parameters 
in  linear  models  under  censorship. 

22.  Sequential  estimation  of  the  regression  parameter  in  the  Cox 
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b)  List  of  Invited  Colloquia,  Conferences  and  Seminars; 


Joseph  C.  Gardiner,  Principal  Investigator 


The  results  of  our  research  efforts  performed  under  Contract 
N00014-17-C-0522  were  presented  at  the  following  professional  gatherings. 

Ic  Symposium  on  Adaptive  Statistical  Procedures  and  Related  Topics, 
Brookhaven  National  Laboratories,  New  York  (1985) 

2.  Institute  of  Mathematical  Statistics,  Eastern  Regional  Meeting, 
Stony  Brook,  New  York  (1985). 

3.  Department  of  Biostatistics,  Columbia  University,  New  York  (1985). 

4.  Department  of  Statistics,  Colorado  State  University,  Fort  Collins 
(1985). 

5.  GTE  Laboratories.  Waltham,  Massachusetts  (1984). 

6.  Department  of  Applied  Mathematics  and  Statistics,  State  University 
of  New  York,  Stony  Brook,  New  York  (1984). 

7.  Statistics  Center,  Cornell  University,  Ithaca,  New  York  (1984). 

8.  Department  of  Statistics  and  Biostatistics,  University  of  North 
Carolina,  Chapel  Hill,  North  Carolina  (1984). 

9.  Department  of  Statistics,  Oklahoma  State  University,  Stillwater, 
Oklahoma  (1983). 

10.  American  Statistical  Association  and  Institute  of  Mathematical 
Statistics,  Joint  Meeting,  Toronto,  Canada  (1983). 

11.  Department  of  Mathematical  Sciences,  State  University  of  New 
York,  Binghamton,  New  York  (1982). 

12.  Department  of  Mathematical  Statistics,  Columbia  University, 

New  York  (1982). 


13.  Biometric  Society  (ENAR),  San  Antonio,  Texas  (1982). 

14.  Special  Topics  Conference  on  Survival  Analysis,  The  Ohio  State 
University,  Columbus,  Ohio  (1981). 

15.  Conference  on  Reliability  and  Maintainability  Office  of  Naval 
Research  and  Army  Research  Office,  Washington,  D.C.  (1981). 
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c)  Description  and  significance  of  research  accomplished  under 
Contract  N00014-17-C-0522 


Joseph  C.  Gardiner,  Principal  Investigator. 

We  shall  briefly  review  here  the  results  obtained  in  the  16 publ ications 
listed  in  (a)  and  indicate  their  contribution  to  statistical  theory  and 
methodology.  Our  discussion  is  presented  under  the  following  rubrics, 
i)  Introduction 

ii)  Sequential  Inference  Procedures  with  ;:ensorship  data 
iii}  Likelihood  ratio  processes  in  progressively  censored  data 
iv)  Biometric  functional  estimation  in  multiple  decrement  models 

i )  Introduction. 

Our  research  effort  has  been  in  the  broad  area  of  statistical  analyses 
with  censored  or  incomplete  data.  The  basic  statistical  model  entails 
observation  of  a  lifetime  X  or  its  competing  censoring  time  Y  ,  which¬ 
ever  comes  first,  for  each  unit  in  a  sample  of  specimens.  The  term 
"lifetime"  is  given  a  broad  interpretation  and  includes  such  examples  as 
the  time  to  failure  of  a  mechanical  or  an  electronic  component  as  in  a 
reliability  study,  or  the  survival  time  of  a  laboratory  specimen  as  in 
a  clinical  study  or  the  duration  to  some  specific  end  point  as  in  a 
demographic  or  epidemological  investigation.  In  these  situations  it  is 
often  the  case  that  the  duration  variable  X  of  interest  is  detered  from 
complete  observation  due  to  the  presence  of  censorship.  Examples  include 
withdrawals  from  study  in  clinical  and  iatrogenic  trials  and  the  partial 
removal  of  units  from  testing  in  a  reliability  study.  The  censoring 
time  Y  may  be  stochastic  or  fixed,  but  the  investigator  has  to  be 
content  to  work  with  the  data  Z=min  (X,Y)  and  the  identifier  6=1  or  0 


according  as  whether  failure  preceded  censoring  or  not.  For  a  sample  of 


n  specimens  therefore,  we  have  the  data  l<i<n}.  In  the  random 

censorship  model  one  assumes  the  X^,Y^  are  independent  random  variables, 
that  is,  the  censoring  mechanism  is  noninformative. 

In  several  longitudinal  studies  the  vat^iables  {Z^. ;  lfi<n}  are  them¬ 
selves  time-ordered;  the  first  observation  to  be  recorded  is  the  smallest 

★ 

one  Z^j  together  with  its  identifier  which  is  1  or  0  according 

as  Z^j  is  an  uncensored  or  censored  observation.  This  is  followed  by 

the  next  smallest  observation  Z^2)  {together  with  62)  and  so  on  until 
* 

the  pair  is  recorded  last.  It  is  often  the  case  that  prolonged 

observation  may  be  necessary  before  the  entire  sample  of  observations 
★ 

{(Z(i),6n):  has  been  recorded.  Limitations  on  time  or  economic 

considerations,  let  alone  ethical  reasons,  preclude  extended  periods  of 

experimentation  and  therefore  curtailment  of  observation  at  an  intermediate 

stage  is  generally  warranted.  If  observation  is  ceased  at  the  k^^-th  stage 

where  k^e{ 1 ,2, . . .  ,n}  may  be  a  stopping  time--an  integer-valued  stochastic 

variable  such  that  [kp=t]  depends  only  on  information  l<i<t} 

gathered  up  to  that  stage--then  the  investigator  has  to  deal  with  the 

★ 

recorded  data  l<i<k^}.  Such  data  will  be  called  progressively 

censored  data  and  statistical  analyses  based  on  them  are  necessarily  more 
complicated  and  subtle  due  to  ordering  aro  dependence  among  the  variables. 

The  basic  statistical  model  we  have  described  here  can  be  extended  to 
a  more  general  competing  risks  or  multiple  decrement  model  where  data  of 
the  type  (X,J)  are  gathered.  Here  once  again  X  denotes  the  duration 
variable  to  some  end  point,  say  failure,  ■.-.'■Me  J  labels  the  cause  of 
failure  or  decrement.  The  random  censor£’'ip  model  may  be  regarded  as  a 
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speclal  case  in  which  there  are  two  causes  of  decrement,  true  failure  or 
censoring,  and  X  is  replaced  by  Z  as  before.  In  the  terminology  of  a 
competing  risks  model  the  variables  X,Y  are  called  latent  variables. 

i i )  Sequential  Inference  Procedures  with  censored  data. 

The  consideration  of  large  sample  si7es  is  often  inappropriate  in 
several  reliability  and  lifetesting  situations  where  economy  of  time  and 
high  per  unit  costs  of  units  preclude  implementation  of  statistical  pro¬ 
cedures  which  require  relatively  large  sample  sizes  for  their  proper 
utilization.  In  this  respect  our  publications  (1,2,4,17)  address  the 
question  of  efficient  estimation,  in  a  sense  to  be  made  precise  here  later, 
of  some  functionals  of  the  underlyiriy  failure  distribution.  To  describe 
our  results  suppose  we  have  placed  n(>l;  homogeneous  units  on  test  and 
record  as  time  elapses  the  true  failure  times  X^  or  the  censoring  time 
Y.  ,  whichever  comes  first  for  each  unit.  We  wish  to  estimate  some 

functional  e=0(F)  of  the  failure  distribution  F(*P[X>‘])  based  on  the 

* 

accumulated  progressively  censored  data  stage 

,  where  k^  is  an  appropriate  stopping  time.  Let  b(>0)  be  the  per 
unit  cost  of  an  unit  placed  on  test  and  c(>0)  be  the  per  unit  cost  of 
time  or  test.  At  stage  k  ,  the  total  time  on  test  (TTT)  V  .  is  given 

~  n  j  In 

by 

’'n.r  • 

Suppose  n=n(c)  units  are  placed  on  test  and  9  .  is  the  current  estimate 

11  y  In 

2 

of  0  while  the  estimation  error  is  measured  by  squared  error  (0  |,-6)  . 

n )  K 

Then  we  may  be  taking  the  loss  incurred  in  estimatinq  9  by  0  .  at 

n )  K 

stage  k  as 

<2-2)  ''n.k 

where  a(>0)  is  a  known  tonstant.  Then  the  natural  objective  is  to 
minimize  the  expected  loss  (  risk),  with  respect  to  k,  for 


given  (n,b,c).  We  have  investigated  this  problem  in  (4)  where  we  assume 
e  is  the  exponential  mean  failure  time  and  the  censoring  distribution  is 


unknown  but  independent  of  the  failure  distribution  G.  The  maximum  likeli¬ 


hood  estimator  of  e  at  stage  k  is  e 


n  |,/6-  ,  where  z  6.  . 

n,k  n,k  n,k  *  n,k  ^ 


Assuming  lim  cn  =a  e(0,®),  b=pc  it  is  shown  that  along  sequences  {k„} 
c+0  " 

for  which  n“^k^->-X,  Xg(0,1]  ,  the  risk 


(2.3)  .  ={a6^/b  )k"^+pcn+cb  k„+o(k”^)  where  b  =x“V!|  ^^^H(x)dx  ; 

FijK^  An  Ann  A  o 

H'{o)b(o)=-l  and  H=FG  .  In  this  situation  the  optimal  choice  of  is 
the  integer  k°  closest  to  (ae'^/cbn^,  if  a  >a0'’(EZ)"^  and  n  otherwise. 

U  A 

In  the  later  case  observation  continues  until  all  units  have  responded. 

It  can  be  shown  that  an  be  optimized  only  along  sequences  {k^}  for 

which  n'^k  +X,  Xe(0,l],  As  oO  notice  that  n“^k°-»-A  ,  provided  X  is 
1 1 

a  solution  of 


(2.4) 


/q  ^^^^H(x)dx=(aeW)^ 


Since  both  x,e  are  unknown  the  specification  of  the  optimal  stopping 

number  k°  cannot  be  made.  We  explore  an  analogous  sequential  procedure 

which  ceases  observation  at  the  random  stopping  stage  and  consequently 

estimates  e  by  •  Defining 

’  c 

(2.5)  N^=  rmin{k<  n-1:  >  (a/c)V^^,^} 

In  ,  if  no  such  k  exists 

it  is  shown  that  the  sequential  estimation  scheme  (N  ,e  )  behaves  as  good 

as  the  "fixed-sample"  scheme  (k°,  e^o)  with  regard  to  their  risk.  If 

c 

*  o 

w  the  risk  of  the  sequential  procedure  and  Ri  is  the  correspond 
c  n  j  c 

ing  optimal  risk  of  the  procedure  (k°,  e  o)  then  R*/R°->1  as  c+o,  that 

c  He  c  c 

is, our  scheme  is  asymptotically  risk-efficient. 
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Furthermore  N  /k°-*-l  and  e  properly  normalized,  has  the  limiting 
normal  distribution. 

Our  publication  (4)  is  the  first  attempt  in  obtaining  a  risk  efficient 
time  sequential  scheme  for  the  estimation  of  the  mean  exponential  failure 
time  when  the  observations  are  subject  to  censorship. 

The  results  of  Sen  (1981),  in  which  censorship  is  absent,  can  be 
obtained  as  a  special  case.  The  quintessence  of  our  arguments  is  the 
derivation  of  (?.3).  In  the  absence  of  censoring,  §  ^  and  since 

Vn  1^  of  is  then  the  sum  of  k  independent  exponential  random  variables 

with  mean  e,  R  .  takes  a  particularly  simple  form.  The  technical  machinery 

11  j  K 

brought  to  play  in  the  case  of  censorship  involved  certain  inequalities 
for  martingales  both  in  discrete  and  continous  time,  convergence  of  sequences 
of  dependent  random  variables  and  the  moments  of  normalized  empirical  pro¬ 
cesses.  We  believe  that  some  of  these  tools  are  also  new. 

In  publications  (1)  and  (3)  the  parallel  problems  of  risk-efficient 
estimation  of  a  and  fixed  width  interval  estimation  are  analysed.  In 
(1)  the  loss  structure  is  simpler,  n'®^®n  n”®^  focus  on 

minimizing  the  recruitment  sample  size  n=n(b)  of  units  for  a  given  per  unit 
cost  b  .  In  (3)  the  objective  is  to  obtain  a  confidence  interval  for  e 
of  given  width  2d  having  a  prescribed  coverage  probability  (1-a).  Since 

„  is  the  maximum  likelihood  estimate  of  e  based  on 
n,n  n,n  n,n 

®n  n  ’  ^®n  limiting  normal  distribution  N(o,o  )  a  natural 

candidate  for  a  confidence  interval  for  e  of  width  2d  is 
In=[§n  ®n  taking  a  sample  size  n=n(d)  of  at  least 

(Zo^/2a/d)  we  can  guarantee 

lim  PQ[I_Ee]>l-a  for  all  6 
y  n  — 

-►0 
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However,  since  a  is  unknown,  depending  on  e  and  the  unknown  censoring 
distribution,  the  specification  of  n(d)  cannot  be  made.  We  propose  a 
sequential  scheme  in  (3)  that  takes  N=N(d)  observations  such  that  the 


resulting  confidence  interval  I, 


is  asymptotically  consistent  and 


efficient,  that  is  as  d-^o  and  all  e. 


lim  RgCeelj^]  =  1-a 


lim  EQ{N(d)}/n(d)=l  . 

If  the  underlying  failure  distribution  is  also  left  unspecified,  then  problem 
of  efficient  estimation  of  the  mean  failure  time  becomes  considerably  more 
difficult.  We  have  addressed  this  problem  in  (2)  and  in  the  case  of  median 
failure  time  in  (5)  and  (17).  The  estimates  chosen  here  are  functionals 
of  the  product  limit  estimator  of  F. 

The  techniques  developed  here  hold  much  promise  in  solving  some  related 
problems  which  we  shall  examine  in  our  future  research. 

In  addition  to  the  optimality  property  of  asymptotic  risk-efficiency 
in  the  point  estimation  case  and  of  asymptotic  consistency  and  efficiency 
in  the  interval  estimation  case  an  investigation  of  an  higher  order 
of  optimality  becomes  necessary  in  certain  situations.  The  basic  problem 
of  efficient  estimation  of  the  median  failure  time  is  taken  up  in  (16,  18, 
19).  A  bounded  length  interval  estimation  scheme  is  one  in  which 
sampling  is  curtailed  upon  achievement  of  a  desired  prespecified  width 
of  the  random  interval  under  utilization.  The  fixed-width  design  described 
earlier  always  maintains  a  constant  (fixed)  width  throughout  sampling. 


'  w 
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In  publication  (18)  we  find  that  both  these  procedures  are  asymptotically 
consistent  and  efficient.  Therefore  theoretical  development  of  expansions 
for  the  coverage  probability  and  of  expected  sample  size  will  be  needed 
to  differentiate  between  these  schemes.  A  simulation  study  under  investi¬ 
gation  at  present  has  shed  some  light  as  to  which  procedure  will  be 
better.  The  current  simulations  carried  out  by  the  investigator 
point  to  the  bounded  length  procedure,  in  a  particular  setting,  as 
appropriate  for  moderate  sample  sizes.  The  problem  is  being  investigated 
further. 


i i i )  Likelihood  ratio  processes  in  progressively  censored  data. 

The  property  of  local  asymptotic  normality  (LAN)  for  a  family  of 
probability  measures  is  of  fundamental  importance  in  the  asymptotic  theory 
of  estimation  and  hypothesis  testing.  A  comprehensive  examination  of  sets 
of  conditions  ensuring  the  LAN  property  for  various  families  of  probability 
distributions  can  be  found  in  LeCam  (1960).  Numerous  applications  are 
available  in  Roussas  (1971)  and  Ibragimov  and  Has'minskii  (1981).  Attention 
has  been  focussed  on  families  of  distributions  connected  with  a  sequence 
of  independent  random  variables  or  those  associated  with  a  homogeneous 
Markov  chain.  For  the  progressively  censored  data  the 

variables  are  neither  independent  nor  identically  distributed.  In  (9) 


this  LAN  property  is  derived  for  a  wide  class  of  parametric  failure  distri¬ 
butions  and  some  applications  to  estimation  are  described.  Suppose  the 
underlying  variable  Z  has  distribution  involving  the  parameter  e  and 
corresponding  to  form  the  likelihood  ratio  statistics 


where  p(*;e)  denotes  the  joint  density  of  Z^|^j  =  (Z^j,...,Zjl^j)  .  As 

described  earlier  expermentation  may  be  curtailed  at  a  (random)  stage 

depending  on  the  accumulated  data  and  so  it  is  instructive  to  examine  the 

stopped  statistics  A„  ^  .  Under  certain  conditions  of  regularity  on  the 

’ 

underlying  response  distribution  we  obtain  with  0^=6+un,"*, 


(2.6)  (^)(0,e^)=exp{uJ^(e)W(t)-iu^tJ^(e)  + 

where  t+k^(t)  ,  te[0,l]  are  certain  (random)  integer  valued  functions. 
t>W(t)  is  Brownian  Motion  in  the  space  D[0,1]  and  a  remainder 
negligible  in  probability.  The  expansion  (2.6)  may  be  called  the  LAN 
property  and  its  importance  lies  in  the  fact  that  tests  of  hypothesis 
of  the  likelihood  ratio  type  (for  contiguous  hypotheses)  can  be  constructed 
and  their  properties  determined  by  well  known  properties  of  the  process  W. 
In  (9)  we  have  extended  to  the  progressively  censored  case  some  results  in 
Ibragimov  and  Has'minskii  (1972,  75)  Phillipou  and  Roussas  (1973)  and 
Inagaki  and  Ogata  (1977).  For  a  special  class  of  survival  distributions 
weak  convergence  of  progressively  censored  likelihood  ratio  processes  in 
D[o,®]  is  illustrated  in  (10). 


iv)  Biometric  functional  estimation  in  multiple  decrement  models. 

In  many  survival  analyses,  as  in  clinical  trials  and  lifetesting. 


estimation  of  the  underlying  survival  distribution  is  of  fundamental 


importance  since  it  summarizes  the  mortality  experience  of  the  items  or 

specimens  under  study.  In  the  case  of  censored  data  the  estimator  of  the 

failue  curve  t->^F(t)(*P[X>t])  of  common  use  is  the  celebrated  product- 

limit  estimator.  Based  on  the  data  {(Z^. ,6^):  lfi<n}  this  estimator  can 

be  expressed  in  the  form  n  »  n 

n  N(Z.)  6.«1J 

(2.7)  .Jj  h+N(Z.)^  '* 

where  the  exponent  denotes  the  indicator  function  of  the  event  in  parenthesi 
n 

Here  N(t)  =  I  [Z.  >  t]  is  the  number  of  items  at  risk  at  time  t.  The 
i  =  l  ^ 

properties  of  t->-F^i(t)  has  been  the  subject  of  intensive  research  in 
recent  times  (Gill  (1983)).  In  (8)  we  have  addressed  the  qi  jstion  of 
estimation  of  the  failure  curve  when  the  investigator  has  at  his  disposal 
only  progressively  censored  data  of  the  type  ^(^(■i)» 
approach  adopted  here  is  a  Bayesian  one  where,  by  imposing  a  Dirchlet  pro¬ 
cess  prior  on  F  and  using  a  weighted  squared-error  loss,  an  appropriate 
generalization  to  the  progressively  censored  case  of  an  estimator  of  F 
given  by  Susarla  and  van  Ryzin  (1976)  is  obtained.  We  also  note  that  our 
estimator  in  (8)  encompasses  both  ^ixed  and  random  censorship.  It  includes, 
of  course,  the  cases  in  which  the  complete  sample  is  observed  (kj^=n)  , 
an  extension  of  an  estimator  of  Ferguson  (1973)  when  no  censoring  is  pre¬ 
sent  and  the  product-limit  estimator  (2.7).  The  asymptotic  theory  of  our 
non-parametric  Bayesian  estimate  of  F  is  discussed  in  (7). 

As  noted  in  the  Introduction  the  competing  risks  model  entails  obser¬ 
vation  on  data  of  the  type  (X,J)  where  X  plays  the  role  of  a  duration 
variable,  such  as  failure  time  and  0  an  index  labelling  the  cause  of 
failure.  Central  to  some  epidemological  and  demographic  studies  is  the 


6- 


evaluation  of  certain  incidence  rates  (fertility  rates,  mortality  rates, 
etc.)  which  may  be  expressed  in  terms  of  the  failure  distribution 
t>F(t)=P[X>tj  and  the  cause-specific  failure  rate  functions  {g^}  given  by 

(2.8)  g^(t)  =  lim  P[t<X<t+At,  J=i|X>t]/At. 

At-*-o 

Thus  g^  may  be  interpreted  as  the  instantaneous  probability  of  failure 
by  decrement  "i"  at  time  t,  given  survival  up  to  t  from  all  risks.  On 
the  assumption  that  c(>l)  mutually  exclusive  risks  are  operative  it  is 
easily  seen  that 

(2.9)  F(t)=exp  (-  g.(x)dx). 

i=l  °  ^ 

The  occurrence-exposure  rate  corresponding  the  decrement  "i"  is 
defined  as  the  functions  a^(t)  where 

(2.10)  a.(t)=/J  gi(x)F(x)  dx//jF(x)  dx 

This  can  be  re-expressed  as  PLX<t,  J=i]/E(Xt)  and  therefore  o^(t)  has 
been  interpreted  as  the  mean  mortality  risk  by  decrement  i  over  the 
duration  [0,t]  .  In  (14)  the  estimation  of  the  {a^}  is  addressed  and 
the  asymptotic  distribution  of  empirical  occurence-exposure  rates  is 
obtained.  It  is  noted  that  in  (14)  our  formulation  is  only  in  terms  of 
observable  entities.  We  do  not  introduce  conceptual  failure  times,  one 
for  each  type  of  decrement,  and  therefore  no  ad  hoc  assumptions  are 
necessary  on  their  joint  distribution.  This  averts  the  identifiabi 1 ity 
problems  that  arise  in  this  connection.  An  interesting  consequence  of 
our  theory  is  that  the  empirical  occurrence-exposure  rates  are  in  general 
asymptotically  correlated  with  the  sign  of  this  correlation  related  to 
the  monotonicity  of  the  cause-specific  failure  rates  of  the  individual 
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risks  (when  all  risks  are  operative).  We  have  cited  examples  where  this 
correlation  is  zero,  thereby  ensuring  at  least  asymptotically,  independence 
between  rates.  In  (15)  a  more  elaborate  view  is  taken  by  examining  the 
behavior  of  the  process  of  empirical  occurrence-exposure  rates  in  the 
space  D^[o,«). 
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Concluding  Remarks 

The  principal  investigator  expresses  his  thanks  to  the  Office  of 
Naval  Research  for  the  support  and  interest  extended  in  the  research 
performed  under  this  contract. 

The  statistical  methodologies  developed  under  this  contract  have 
provided  practical  and  useful  extensions  of  the  literature  on  the 
analyses  of  failure  time  data.  Furthermore,  this  has  led  to  new  results 
in  the  treatment  of  censored  sequential  designs  in  survival  analyses. 

Our  current  program  of  research  will  address  various  important  early 
decision  rules  including  curtailed  testing  procedures  that  will  have 
impact  on  statistical  practice  in  several  longitudinal  studies. 
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